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Abstract 

The adsorption of charged colloids (macroions) onto an oppositely charged planar substrate is 
investigated theoretically. Taking properly into account the finite size of the macroions , unusual 
behaviors are reported. It is found that the role of the colons (the little salt-ions carrying the 
same sign of charge as that of the substrate) is crucial to understand the mechanisms involved in 
the process of macroion adsorption. In particular, the colons can accumulate near the substrate's 
surface and lead to a counter-intuitive surface charge amplification. 
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I. INTRODUCTION 



Whereas the bulk behavior of homogeneous (chargedi*^ and uncharged^) colloidal sus- 
pensions is rather well understood, the situation for its inhomogeneous counterpart, such as 
that emerging in an adsorption process, is less clear. Potential applications of adsorption of 
charged colloidal particles (macroions) can vary from technological processes such as surface 
coating^ to biological material problems. - From a fundamental point of view, the tremendous 
long-ranged Coulomb interaction that sets in represents a formidable theoretical challenge. 
Consequently, a deeper understanding of the phenomenon of macroion adsorption is justified 
and needed. 

On one hand, experiments^*^ and the well known mean field Gouy-Chapman theory^i^ 
seem to nicely agree for the ion distribution of an aqueous monovalent electrolyte near pla- 
nar charged interfaces, as long as non-specific forces as well as excluded volume effects are 
negligible. On the other hand, if the solution contains highly multivalent and/or large sized 
ions, then the Gouy-Chapman theory may severely qualitatively fail. A crucial missing in- 
gredient in this theory is the inclusion of the finite size of the macroions that can lead to 
non-trivial phenomena, such as substrate's surface charge reversal, already with monovalent 
ions.— liiii^ii^ The lateral macroion-macroion electrostatic correlations that are also absent 
in a mean field theory can be attenuated (and even become marginal) at sufficiently high 
salt content, in contrast to excluded volume effects. In the past, Gonzales- Mozuelos and 
Medina-Noyolar^ used integral equations to address the problem of macroions near repul- 
sive/attractive charged walls interacting via an effective Yukawa potential. More recently, 
Neta^ investigated thoroughly and analytically the behavior of macroions near charged in- 
terfaces, based on the Debye-Hiickel (DH) theory, but ignoring its finite size. Thereby, 
the striking effect of surface charge amplification advocated here could not be captured 
by those approaches.— It is only very recently, that this phenomenon was reported by 
Lozada-Cassou and coworkers.—*^ For a size-asymmetrical electrolyte,^^ surface charge am- 
plification was identified by solving numerically the modified Gouy-Chapman (where cations 
and anions have different distances of closest approach to the interface). Intimately related 
to our work, the phenomenon of surface charge amplification in presence of macroions was 
also found by applying a sophisticated hypernetted chain/mean spherical approximation 
(HNC/MSA) integral equation.i^ 
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In this work, we present a simple model, where the finite size of the macroions is taken into 
account, to reveal the mechanisms governing macro ion adsorption. Although our approach 
is less accurate than the one used by Lozada-Cassou and coworkers,— it presents the nice 
advantage to be analytical and very intuitive. The basic driving force of surface charge 
amplification is that (spherical) macroions tend to be surrounded by its counterions over its 
whole surface in a uniform manner. This corresponds actually to the old classical Thomson's 
sphere problem.— As long as the strength of the surface charge density of the oppositely 
charged substrate is low enough, a finite number of counterions of the macroions should stay 
in the vicinity of the interface (see Fig. [1]), leading to a surface charge amplification. Our 
paper is organized as follows: Our (modified) DH theory is explained in Sec. [ITI The results 
are presented in Sec. Illll and followed by a brief summary (see Sec. lIVp . 

II. DEBYE HUCKEL THEORY 

A. Linearized Poisson-Boltzmann equation 

Our electrostatic model (see Fig. [1]) resembles that employed by Spitzer- who studied 
monovalent size-asymmetrical ions near a wall. The negative substrate's surface charge 
density is denoted by ctq. The macroions carry a positive central charge +Zme with e 
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z- distance from charged surface 



FIG. 1: Schematic view of our model setup. The macroions are characterized by a distance of 
closest approach z = a to the charged surface, leading to two screening strengths kq and k for 
< z < a and z > a, respectively. 
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representing the usual elementary charge and Zm its valency. Excluded volume effects are 
taken into account via the distance of closest approach a to the charged interface (see Fig. [1]). 
Our model corresponds somehow to the minimal correction for macro ion' size effect. Having 
to deal with a finite concentration of macroions at contact with a reservoir, whose bulk value 
is given by n^,, electroneutrality requires a bulk macroion's counterion concentration Z^n^ 
(assuming monovalent point-like anions, i.e. Z_ = — 1). Those anions will be referred to as 
the coions of the substrate. Additional (point-like) monovalent counterions (Z+ = +1) and 
colons are also considered with a bulk salt concentration Ug. 

An intuitive and widely used way to connect self-consistently the electrostatic potential, 
"^{z), to the total charge density, p{z) ^ J2a''^aZaeexp{—pZae'^) (with a standing for the 
ionic species), is provided by the so-called Poisson-Boltzmann equation that reads: 

A^(z) = ^ [-2n3sinh(e/3^) - Z„ra„exp(e/3^) 

6Q6r 

+Zmnm exp{-ZmeP'^)Q{z - a)] , (1) 

where (^r) is the vacuum (relative) permittivity, /? = l/iksT) is the reduced inverse 
temperature, and is the usual step (or Heaviside) function. This non-linear differential 
equation ([T]) can only be solved numerically. However a linearization of Eq. ([T]) [i.e., DH 
approximation], valid for \e(3'^\ <^ 1 when < z < a and for \Zme/3'$\ <^ 1 when z > a, 
permits an analytical treatment that is going to be discussed. 

For the first diffuse region (0 < 2; < a) made up uniquely of the little monovalent ions 
(see also Fig. [1]), the DH equation is given by Au{z) = (/t^ + k^)u + with = Snisns 
and = ATriBZmnm, where we have introduced the Bjerrum length £3 = ^T^eS UbT ^"^^ 
dimensionless variable u = e/3\l/. The corresponding solution reads: 

u{z) = de-^"' + C2e^''' - ^ {0 < z < a), (2) 

where Kq = + = 47r£B(2ns + ZmXim)- Ci and C2 are integration constants that are 
going to be determined after having applied the suitable boundary conditions. 

For the second diffuse region {z > a) containing all the ions (including the macroions - 
see also Fig. [T]) the DH equation reads Au{z) = k^u with = Kq + Z^i^I = inisl'^ns + 
Zm[Zm + l)'^m]- The physically sound solution with vanishing electric field at 2; — +cx3 
corresponds to: 

u{z) = Cse"'"' {z>a), (3) 
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where C3 is a third integration constant. 



B. Boundary conditions 

Our electrostatic model system is completely characterized once the integration constants 
Ci, C2 and C3 appearing in Eqs. ([2]) and ([3]) are specified. To do so, we apply the three 
following boundary and/or matching conditions: (i) The Gauss' law applied at the charged 
interface 2; = requires u'(0) = |, where h = ^^^^|^^ is the so-called Gouy-Chapman 
length, (ii) The zero intrinsic surface charge at z = a imposes the continuity of the electric 
displacement D = e^SrE and hence also that of the electric field E [i.e., u'{z a~) = 
u'{z a^)]}^ (iii) The continuity of the electrostatic potential dX z = a requires u{z — >■ 
a^) = u{z —>■ a+). The resulting set of three equations can be readily solved and yields: 
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Note that in the limit of point-like macroions (koO 0) one recovers the well-known result 
C3 — ^ — At this stage all the relevant observables of the system can be in principle 
obtained within the framework of the DH theory. 



III. RESULTS AND DISCUSSION 



A pertinent quantity that is appropriate to characterize the strength of the macroion 
adsorption is provided by the contact potential of interaction Um- The latter corresponds 
to the external work accomplished upon bringing a macroion from infinity {z = +00) to 
contact {z = a): Um = — J!^^ ZmeE^dz (with E^ = — ^ denoting the z-component of the 
electric field). With the help of Eq. we get = Zmu{a) = Z^CaC"'^". Using then Eq. 
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for C3, we obtain the following expression for Um- 



PUr, 



2 

T7 — T + ^ tanh (Kna) 



/to6 cosh (Kga) 



— + tanh [Koa) 



(7) 



As can be seen from Eq. ([7]), Um depends on many parameters, such as Zm, b, a, k and Re- 
making its understanding rather difficult. Nonetheless, Eqs. ([2]) - ([7]) suggest that more 
insight on the adsorption behavior can be gained by considering two limits: (i) the high 
screening regime {KQa ^ 1) and (ii) the weak screening regime (koa ^ 1). Without loss of 
generality, we are going to explore these two limits. 



A. High screening regime 

In the high screening regime {nQa ^ 1), Eq. ([7]) can be approximated by: 

f3Um ^ l~ (8) 

where the relation = ^(k^ — Kq) has been used. Interestingly, in this regime, Um depends 
only on the screening strength contrast k/kq- One can also conveniently express U as a 
function of the bulk concentrations rtm and ra^, which reads 



(3U ~ 1 j'^^s ~\~ Zm{Zm -\- '^)nm 

Equation iQ shows that Um is even independent of the Bjerrum length is, which means 
that the effective wall-macroion attraction is entropically (or depletion) driven. To better 
understand this phenomenon, we have sketched on Fig. [2] the profile of u{z) as well as 
profile of the reduced (dimensionless) electric field E*[z) = —^u'{z) [with E*{z = 0) = —1] 
for a set of typical parameters of charged colloidal suspensions: a = 0.5/im, b = 10~^/xm 
(corresponding to |cro| ~ 0.036Cm~^ for an aqueous solvent), Zm = 10"^, Ug = 10~^M, a 
(fictive) equivalent volume fraction (pm = Um^na^ set to 10^^, (i.e., Um ~ 3.1714 x 10~^^M 
with M standing for mole per liter units). Thereby we have K,Qa ^ 16.4. The profile of E*{z) 
reveals an unusual non-monotonic behavior near contact. This is a direct consequence of an 
accumulation of excess" colons in the macroion depleted zone {z < a) leading to a weaker 
screening. 

In order to further characterize the adsorption behavior in the high screening regime, we 
examine Um as predicted by Eq. ([7]) and Eq. ([9]) for some typical experimental values of 
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FIG. 2: Profiles of the reduced electrostatic potential u{z) (solid lines) and the reduced electric 
field E*{z) = — |m'(z) (dotted lines) around contact {z = a = 0.5^m) calculated from Eqs. ([2|) and 
([3]). The inset corresponds to a wider z-range. The following parameters were chosen: b = 10~^/xm, 
is = 7.14 X 10"^/im, = 10^, n, = IQ-^M, and = lO'^. 
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FIG. 3: Reduced contact potential of interaction as a function of the molar bulk salt 

concentration Ug in the strong screening regime with a = 0.5/im and (pm = W^^. From left to 
right: Zm = 10^, 10^ and 10^. The solid lines were generated using Eq. Q. The dotted lines were 
generated using Eq. ([7|) with b = Inm and £b = 0.714nm. For Z,„ = 10^ the two curves are in 
practice identical. 

the parameters k and kq of charged colloidal systems. In Fig. [3l —pUm = P\Um\ is plotted 
against the bulk salt concentration for different prescribed values of = 10^, 10^ and 
10^. Only values of /9|t/m| smaller than unity are shown such as to explore its behavior 
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where the DH approximation is vahd. Figure [3] shows that beyond a certain threshold of 
sah content that is Z^-dependent, no adsorption occurs as signaled by a nearly zero value 
of Urn- At relatively low enough salt concentration, adsorption is favored where Um increases 
with growing Z^- The latter point can be clearly understood by noticing that for typical 
values of the parameters of charged colloidal systems iZ^"^ <^ 1 and Z^ ^ 1) we have 
f3Um ^ 1 - ^1 + Z^^. Moreover, in the limit Z^^ < 1, the simple following relation 
/3f/„ ^ -Zl^ holds. 



B. Weak screening regime 

We now address the weak screening regime characterized by hqq <^ 1. In this situation 
Eq. ([7]) becomes: 



2Z 



m 



Kb 
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kq \ nb Ik 



(10) 



where the dispersion relation k? = + Z^i^^c been used. Equation (ITU]) reveals that Um 
varies affinely with macroion size a. Depending on the sign of ( ^ — 1 ) — 1, Um may 
either decrease or increase with a. For ^^y- ^ 1 we obtain the following limit behavior: 



k2 



lim (3Um ^ - [I - ^] f^a. (11) 

In this "entropic" limit \Um\ increases linearly with a, meaning that a finite wall-macroion 
attraction persists even at vanishing surface charge density ctq due to excluded volume effect 
(a 7^ 0). Noticing that 

Kb 

E*{a) = -—PUm, (12) 

we deduce that the strength of the electric field at contact |i?*(a)| becomes larger than unity 
when 2^ is sufficiently large [especially true for Eq. f|TT]) ]. This feature corresponds to an 
electric field (or surface charge) amplification at contact. 

Let us consider some typical experimental values of parameters representative of charged 
micellar systems (see caption of Fig. Hj). Thereby we have K^a ~ 0.177. To access the 
mechanisms of wall-macroion attraction in the weak screening regime, we have plotted the 
profiles of u{z) and E*{z) in Fig. H] for various values of b. Strikingly, in the macroion 
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FIG. 4: Profiles of the reduced electrostatic potential u{z) (solid lines - from top to bottom: 
b = 500, 300, 50nm) and the reduced electric field E*{z) = — (dotted lines - from top to 
bottom: b = 50, 300, 500nm). The inset corresponds to a magnification of E*{z). The following 
parameters were chosen: a = 5nm, = 0.714nm, = 10, Ug = lO^^M, and (j)^ = 10^'^. 

depleted zone, the electric field gets monotonically amplified (from z = up to contact z = 
a = 5nm) at (very) large b [corresponding to poorly charged interfaces (|cro| < 10~^Cm~^) 
- see Fig. H]. This phenomenon is again due to an accumulation of "excess" colons in the 
macroion depleted zone, that leads here to a net surface charge that is more negative than 
ctq. It is important to remark, that surface charge amplification does not necessarily involve 
a subsequent charge reversal, as clearly indicated in Fig. |H In fact this charge-amplification 
phenomenon was already observed, although uncommented, in computer simulations [see 
Fig. 4(b) in Ref.^ and Fig. 3 in Ref.-Si]. 

IV. CONCLUDING REMARKS 

To summarize, we have studied analytically within the framework of the Debye-Hiickel 
theory the (weak) electrostatic adsorption of macroions at oppositely charged planar sur- 
faces. Taking into account the crucial role of the finite size of the macroions, our model 
reveals non trivial adsorption driving forces. In the strong screening regime, the wall- 
macroion attraction strength at contact is exclusively governed by the screening strength 
contrast k/kq. In the weak screening regime, the wall-macroion attraction strength can 
either decrease or increase with macroion size, depending on the surface charge density of 
the substrate. In particular, at sufficiently small surface charge densities an effective electric 
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filed (or equivalently an effective surface charge) amplification sets in. All these adsorp- 
tion mechanisms have a common feature, namely, the accumulation of excess colons in the 
macroion depleted zone. The latter mechanism also explains the surface charge amplification 
recently reported by molecular dynamics simulations for DNA-dendrimer complexation.— 
Our findings could be experimentally verified by employing an (extra) ultra fine dispersion 
of charged particle tracers (with a smaller size than the macroions and such as to nearly not 
modify the screening strengths /t and kq), whose density profile should reveal the electric 
filed' one in the macroion depleted zone. 
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